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Abstract. Observations of X-ray cavities formed by powerful jets from AGN in galaxy cluster cores 
are commonly used to estimate the mechanical luminosity of these sources. We test the reliability of 
observationally measuring this power with synthetic X-ray observations of 3-D MHD simulations 
of jets in a galaxy cluster environment. We address the role that factors such as jet intermittency 
and orientation of the jets on the sky have on the reliability of observational measurements of cavity 
enthalpy and age. An estimate of the errors in these quantities can be made by directly comparing 
"observationally" derived values with values from the simulations. In our tests, cavity enthalpy, age 
and mechanical luminosity derived from observations are within a factor of two of the simulation 
values. 
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INTRODUCTION 

X-ray images of giant cavities in galaxy clusters associated with powerful jets from 
central active galactic nuclei (AGN) suggest that AGN may play an important role in 
the energetics of galaxy clusters (e.g. [10, 1]). The minimum energy required to produce 
cavities, the enthalpy, is often in the range of 10^^ to 10^° ergs [[l|]. The energy required 
to suppress the radiant cooling of the intracluster medium (ICM) below 2 keV is on the 
same order as energy in these cavities [3]. One popular hypothesis that has emerged to 
solve this "cooling flow problem" is that AGNs release sufficient energy into the ICM 
to quench cooling. A common method for estimating the energy released by an AGN 
requires measuring the cavity enthalpy and age [1]. We have produced synthetic X-ray 
observations from a pair of 3D magnetohydrodynamic (MHD) simulations of jets in 
realistic cluster environments. The two simulations tested the observational implications 
of jet intermittency on cavities. 



SIMULATIONS 

We present two simulations of jets in galaxy clusters representing two models of jet 
intermittency [5]. The simulations were computed using a 3D MHD total variation 
diminishing (TVD) code that employs flux constrained transport to maintain divergence- 
free magnetic fields ([tzilsl]). The calculations were done on a Cartesian grid 600 zones 
in X and 480 zones in y and z with Ax=Ay=Az = 1 kpc. Passive relativistic electrons were 
included following the Coarse-Grained finite Momentum Volume scheme 
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FIGURE 1. Synthetic X-ray observations of RE {left) at 158 Myr and 113 {right) at 171 Myr with 
the grids rotated to an incHnation of 45" from 1.5 to 2.5 keV divided by a best fit double-j3 profile to 
emphasize the cavities. Note the presence of ripples in the 113 image. These are produced as the jets are 
toggled on and off. The 113 cavities have a wider aspect ratio as compared to RE. 



The calculations were both initialized with an ambient environment intended to mimic 
many characteristics of observed galaxy clusters. A static gravitational potential was 
included defined by an NFW dark matter distribution [4] with a virial mass of 4.2 
X IO^'^Mq for a virial radius of 2 Mpc and core radius of 400 kpc. The ICM was 
set in hydrostatic equilibrium with a core pressure of 2.4 x 10^^^ dyne cm^^ and 
core density of 8.33 x 10^^^ g cm^-^. A Kolmogorov spectrum of density fluctuations 
with an amplitude of ~10% was added on top of the initial density profile. The ICM 
magnetic field was ~10 /iG and tangled on scales of ~10 kpc. A pair of bi-directed 
jets originated from the center of the computational domain. The jets were launched in 
pressure equilibrium with the surrounding ICM at Mach 30 with respect to the ICM 
sound speed, vjet ~0.1 c. The jet diameter was 6 kpc . Jet plasma was injected onto the 
grid with a density contrast of Picm/ Pjet = 100 compared to the core ICM density. The 
jets contained a ~10 /iG toroidal magnetic field. Two models of jet intermittency are 
tested; a relic model (RE) of jets turned on for 26 Myr and then permanently turned off, 
and an intermittent model (113) of jets toggled on and off every 13 Myr. Jet power in 
both cases was intended to represent typical FRII radio jets. 

Synthetic X-ray observations of both simulations were calculated over a band from 
0.75 to 10 keV as a post processing step using a ray casting technique. Both thermal 
bremsstrahlung emission and inverse Compton scattering of CMB photons are included 
in the emission model. The computational domains were set at a luminosity distance of 
240 Mpc to produce images with a resolution comparable to Chandra. The grids were 
rotated to several inclinations to test for potential projection effects on the analysis. 

PV WORK & MECHANICAL LUMINOSITY 

"Observational" measurements of the minimum energy to produce the X-ray cavities 
visible in the synthetic X-ray observations were made following a procedure similar to 
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FIGURE 2. "Observed" enthalpy, Ho^,, for RE {left) and 113 (right) measured at various times through 
the simulation. Values are shown from observations made with the grid rotated to three inclinations. The 
enthalpy measured directly from the simulation, Hs,„,, and the total energy added to the grid from jets, 
AEfof, are shown for comparison. 



TABLE 1. Several cavity age estimates, enthalpy and mechanical luminosity measured from "observa- 
tion" for both RE and 113 at 158 Myr and 171 Myr respectively. The jet luminosity in the simulation, 
Lsim, is given for comparison. 
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that of Wise et al. and Rafferty et al. X-ray cavities were fit by eye with ellipses 
to define their projected extent and volumes. We then found best fit double-/3 brightness 
profiles to X-ray images excluding the cavity regions. Assuming an isothermal ICM, 
these profiles were inverted to produce an ICM density profile as given in Xue & Wu 
UM- With the assumption that the cavities were inflated in pressure equilibrium with the 
ICM, an estimate of the minimum energy required to produce the cavities, the enthalpy, 
is given hy H = jpV / (7 — 1). Here 7 is 5/3 to match the simulation physics. 

Three commonly used cavity age estimates were calculated from the observations; 
buoyant rise time t^uoy, refill time tr, sound crossing time [3]. The mechanical lumi- 
nosity of the jets was estimated from the enthalpy and an average of these three cavity 
ages as L,nech = H / {t) . 



DISCUSSION 



The measurement of cavity enthalpy from the synthetic observations produced reason- 
able estimates compared to the actual enthalpy in the cavities in the simulation. "Ob- 
served" values tended to be larger than the simulation values by less than a factor of two 
for both RE and 113 as shown in Figure [21 Projection might have had a large effect on the 
accuracy of the measurements. However, two factors counteract one another to produce 
only a small net variation in the measured values; cavities appear smaller and to be in a 
higher pressure environment closer to cluster core at small inclinations. There was little 
or no dependence on the age of the cavities on the accuracy of the measured enthalpy. 
This was largely due to the fact that cavities closer to cluster core, such as young cavi- 
ties, have a higher contrast ratio (21] making it easier to define their extent. The portion 
of older cavities in this regime contain the bulk of the total energy. Cavity ages from 
thu„y and tc varied significantly due to the variation in the projected distance from cluster 
center. The ages estimated from ty have little dependence on orientation as the estimate 
relies only on a measurement of the gravitational acceleration and cavity radius. L„jgc/i 
derived from the enthalpy and an average of the age values was within a factor of two 
of the real jet luminosity across models of intermittency and jet orientation as shown in 
Table [B The energy required to offset cooling in clusters can be characterized as T^pV 
[111]. An approximate factor of two span in r\ in our tests is due largely to uncertainties in 
the measurement of pV. 
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